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We study the dynamics of a self-gravitating scalar field solitonic object (boson star) in the Jordan- 
Brans-Dicke (BD) theory of gravity. We show dynamical processes of this system such as (i) black 
hole formation of perturbed equilibrium configuration on an unstable branch; (ii) migration of 
perturbed equilibrium configuration from the unstable branch to stable branch; (iii) transition from 
excited state to a ground state. We find that the dynamical behavior of boson stars in BD theory 
is quite similar to that in general relativity (GR), with comparable scalar wave emission. We also 
demonstrate the formation of a stable boson star from a Gaussian scalar field packet with flat 
gravitational scalar field initial data. This suggests that boson stars can be formed in the BD 
theory in much the same way as in GR. 
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I. INTRODUCTION 

Gravitational solitonic objects are quite an interesting 
topic in general relativity (GR). A boson star consists of 
massive complex scalar field, and was first discussed by 
Kaup Q and then by Ruffini and Bonazzola [g (for a 
thorough review, see P,0|). They can form stable con- 
figurations having negative binding energy, as a result of 
a balance of the dispersion due to the classical analog 
of the uncertainity principle and the attractive effects of 
gravity. If we include even a small self-interaction term, 
then their maximum allowed stable mass can be close to 
the order of a solar mass ||. It is also speculated that 
they are a form of dark matter that could have been cre- 
ated during a phase transition in the early universe (see 
Frieman et al M). Although we have still no evidence for 
their astrophysical existence, these systems are a good 
model to learn the nature of a strong gravitational field. 

The stability of boson star has also been studied by 
several authors. Lee and Pang J7J discussed ground state 
stability using linearized perturbation theory, and Seidel 
and Suen p) studied their dynamical behavior by evolv- 
ing field equations numerically. In GR, the ground state 
boson star configurations comprise of a stable branch and 
an unstable branch. Upon perturbations, boson stars on 
the stable branch remain in the stable branch, settling 
down into a configuration with a different mass. In the 
process, it emits scalar radiation with some characteris- 
tic normal mode frequencies. On the other hand, stars 
on the unstable branch do not remain these after per- 
turbations. They either disperse completely, form black 
holes or migrate to the stable branch depending on the 
size of the perturbations. These qualitative features are 
also discussed applying catastrophe theory || . Recently, 
Balakrishna, Seidel and Suen |10| studied dynamical bo- 
son stars with the self-coupling term and excited states. 
They found that excited state boson star equilibrium con- 
figurations have branches similar to that of grand state, 



but all branches are unstable. 

In this paper, we study the dynamical behaviour of 
boson stars in the Jordan-Brans-Dicke (BD) theory of 
gravity pM . In particular, we compare the dynamics 
with those in GR. The BD theory is one of the alternative 
theories of gravity to GR, and the most simple and proto- 
type in all the scalar-tensor theories of gravity (see Will 
|12| for a review). The previous experimental test using 
the delay of radar echo in the solar system shows the 
bound of the BD parameter as lubd > 500 I^Jlif , of 
which the infinite limit agrees with GR. This bound is 
also considered to be limited by the direct observations of 
gravitational waves (see pa] and the references there in). 
[Recently, more strict limitations of the BD parameter 
ujbd are reported JL6]], however the results are model 
dependent and we think lubd > 500 is still the generally 
accepted observational limit.] 

So far, boson stars in scalar-tensor gravity have also 
been discussed by Gundersen and Jensen |r7fl and Torres 
p8[ showed the existence of equilibrium ground state bo- 
son star solutions in the BD coupling, and in the three 
different couplings in the scalar-tensor theories, respec- 
tively. Comer and Shinkai |PJ] showed the existence 
of excited state boson stars in both BD and Damour- 
Nordtvedt's quadratic coupling (attractor) model [£0| in 
the scalar-tensor theories. They also discussed the stabil- 
ities of ground state boson stars using catastrophe theory. 
One of our purposes in this paper is to study the stabil- 
ity of boson stars in BD theory. By evolving the field 
equations numerically for slightly or heavily perturbed 
equilibrium data, we clarify the stability of boson star 
both for ground and excited states. 

Our system includes two scalar fields: bosonic matter 
(complex and massive) scalar field and gravitational (real 
and massless) scalar field (which hereafter we call BD 
field). We expect to be able to study the fundamental 
mechanism of the interactions between these two fields in 
their dynamics. Several previous simulations have shown 



emissions of scalar waves in BD |21j |23J or scalar-tensor 
theory [£4J,£5| from a collapse of dust or star model. Our 
second interest is in how much difference appears in the 
scalar gravitational wave emissions during the dynamical 
boson star system between BD and GR: whether they will 
be enhanced or supressed. 

Throughout this paper, we stand at the point to see 
if they in any way differ from boson stars in a detectible 
manner; that is, we impose ojbd > 500 in most simula- 
tions. 

The outline of this paper is as follows. In Sec. O, 
we will introduce the field equations, and basic outline 
of our numerical techniques. With a view to study this 
problem in scalar-tensor theories in the near future [B2| , 
we describe fundamental equations not only in BD theory 
but in general scalar-tensor theories. 

We des cribe the equilibrium configuration briefly in 
Sec. 



Ill 



We plot the sequences of excited states 
configurations in BD theory and also discuss fractional 
anisotropy measurements in this system, both of which 
have not appeared in the references. 



Sec. IV is devoted for stable boson star configurations. 
In GR, under the effects of finite perturbations, these 
stars react by expanding and contracting with emitting 
scalar radiation at each expansion l|. The star looses 
mass and settles to a lower mass configuration. Each 
expansion (contraction) of the core of the star is accom- 
panied by the contraction (expansion) of the radial met- 
ric. We study whether the expansions and contractions 
of the star set the BD field into oscillations and whether 
this results in any measurable gravitational radiation in 
addition to the scalar radiation. Under infinitesimal per- 
turbations, boson stars in GR start oscillating at their 
fundamental quasinormal mode frequencies that remain 
constant and virtually undamped for large periods of 
time. We also study whether this is seen in BD theory. 

Sec. |v| shows the behavior of unstable boson stars. 
In GR, equilibrium boson stars on the unstable branch, 
when perturbed, begin forming black holes or migrating 
to a new equilibrium configuration on the stable branch. 
The excited state configurations of boson stars are ex- 
pected to be unstable. If they cannot lose enough mass 
and settle to the ground state configuration they collapse 
to black holes. We also show the cascading of a BD boson 
star in an excited state to the ground state. 

In Sec. 



VI 



we look at the formation of boson stars 
in BD theory. The collection of bosonic matter localized 
in a region of space is represented by a Gaussian initial 
boson field. The BD field itself is initially set to zero, 
we discuss whether or how an equilibrium config urat ion 
actually forms in BD theory. Finally, in Sec. VII we 
make some concluding remarks. 

We use the units c = 1 and % = 1. This implies that 
the scalar field mass m is an inverse length (actually, the 
inverse Compton wavelength of the scalar particles) and 
the bare gravitational constant G* has units of length 
squared. 



II. THE FIELD EQUATIONS 

In this section, we present our basic field equations, 
boundary conditions, and numerical techniques to solve 
the system. For our future convenience p3], we show 
the field equations not only in the Brans-Dicke theory 
of gravity (uj = ljbd = const.), but also in the general 
scalar-tensor theory of gravity (u> — ^>(4>))- 



A. Lagrangians and field equations 

The action for our system of scalar-tensor gravity cou- 
pled to a self-interacting, complex scalar field in the phys- 
ical, "Jordan-frame" is 



Sj 




tk-^wi&grdnW^ 



g^d^d^ + V (^ty) 



(2.1) 



The gravitational scalar is <j) and lo((J)) is the Jordan- 
frame coupling of <fi to the matter. The complex scalar tp 
(with its complex conjugate being ijj') has mass m and 
is self-interacting through the potential 



TfP' A 



(2.2) 



The strength of the self-interaction, A, is normally taken 
to be positive. 

There is an alternative representation of the action 
above, the so-called "Einstein-frame." The expression 
is given by the conformal transformation 



~ _ „2a(ip) 
9yM — 6 9fiu 



(2.3) 



where a((p) is the functional transformation from <fi to the 
Einstein-frame gravitational scalar ip, 



0- 1 = G*e 2a(v) , 



(2.4) 



where G* is the effective gravitational constant in the 
Einstein-frame. The relationship between u>(4>) and a((p) 
is obtained from 



a 2 = (2LU + 3)- 1 



where 



da 
aft,) = - . 

The action in the Einstein-frame is thus 
1 



(2.5) 



(2.6) 
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16ttG* 



d*Xy/=g[R-2g' ,v d li <pd v (p] 



(2.7) 



It does not deliver GR exactly because the metric g^ is 
not the true, physical metric that encodes the distance 
between spacetime points. However, the Einstein-frame 
does deliver equations that are similar enough to GR that 
we will use it for our calculations. 

The Einstein-frame stress-energy tensor is 

_I e 2a( V ) (d^ffr,/, + 2e 2a ^V(^V)) 9»u ■ (2.8) 

The gravitational field equations for g^ v and ip are 

G>„ = SirG^Tftu + 2dfj,ipd u ip - d T ipd T ipg tlu (2.9) 

and 

V CT VV = ~4naT , (2.10) 

where T is the trace of the stress-energy tensor. The 
matter field equations are 

BV 

V^VV + 2ad T ^d T ip = 2e 2a{ip) — , (2.11) 

dtp 

dV 
V ff V a tp + 2ad T ipd T ip = 2e 2aM ■— . (2.12) 

dip' 

The coupling function a(ip) is given by choosing a the- 
ory of gravity. In this paper, we only consider the Brans- 
Dicke coupling 



a(ip) = 



f- V>o 



^/2u>bd + 3 ' 



(2.13) 



where the parameter u = lobd is constant, which obser- 
vational constraint is known as lord > 500 [[13 14 1. The 
term ip^ represents the asymptotic value of the gravita- 
tional scalar field. 

Because the potential V (ip ip) ^ s a functional of tp^tp, 
it preserves the global U(l) gauge symmetry (tp — ► e ta tp, 
where a is a constant) present in the theory. This sym- 
metry results in a conserved current, whose explicit form 
in the Jordan-frame is 



2 



JH = l e -2a (V)g ^ f^Q^ _ ^Q v ^ 



(2.14) 



This conserved current leads to a conserved charge, which 
is N p , the number of particles making up the star: 



N p = I d 6 Xy 



-gJ f 



(2.15) 



The spacetime considered here is spherically symmet- 
ric, with the Einstein-frame metric taking the form 



ds 2 = -N 2 (t, r)dt 2 + g 2 (t, r)dr 2 + r 2 [d6 2 + i 



sin 



(2.16) 



where N is the lapse function and r is the circumferential 
radius. We dropped the shift vector, since we use a polar- 
slicing condition [£6| in evolution. 

In this coordinate system, the Jordan-frame ADM 
mass Mj is given by 

G*M, = Una £ (1 - l/g rr ) . (2.17) 

The similar Einstein-frame ADM mass Me is 

G*M E = Hm T - (1 - l/g„) . (2.18) 

r — >oo 2 

However, since g rr — e 2a< - v ^g rr and we set a(^oo) = [to 
be discussed below], then the limits on the right-hand- 
sides are equal and therefore Mj = Me = M. 



B. The equilibrium state equations 

The gravitational scalar, which is real, is assumed also 
to be spherically symmetric and static when we solve the 
equilibrium configuration: 



if = ip(r) . 



(2.19) 



As for the matter scalar field, Friedberg et al |27|] show 
that the minimum energy configurations are those for 
which 



i> 



-ait 



$(r) , 



(2.20) 



where f2 is real and positive and <&( r) is real function. 
Their proof (see the Appendix in |27| ) also goes through 
for scalar-tensor gravity, and so we will take ip to have 
this form. 

We will take advantage of scale-invariances of the field 
equations to redefine some of the fields, parameters, and 
the radial and time coordinates: 



mr — » r, y/AirG*® — > $, mN/fl 
A/47tG*to 2 — > A, Qt/m —> t. 



N, 



(2.21) 



Note that the rescaling changes the asymptotic value of 
N, which is now 



lim N(x) = m/Cl . 

r — >oo 

The field equations, then, become 



(2.22) 



d r d r ip 



1 



2e 4a rg 2 V($)]d r ip 

r 

-g 2 {2ae 2a [-V ff $V CT $ + 2e 2a V r (<J>)]}, (2.23) 



d r d r $ = [-- 



2e ia rg 2 V{$)}d r $ - -^$ 



N 2 



d r (g 2 



„ 2r o„dV{§) 



2 f 9 r 

9 (jp T oo 



aV CT $V>], 



d r {N 2 ) =N 2 (rT 11 + 



,9 2 -l. 



(2.24) 

(2.25) 
(2.26) 



where 



N 2 



Too = —%(d r ip) 



9 



N 2 
+e 2a [^ 2 + 1 -L 
9 



(d r <py + 2N 2 e 2a V(<P)l 



-e 2a [(d r $) 2 + |p2$ 2 - 2g 2 e 2a V{<$>)] 



(2.27) 
(2.28) 



The boundary conditions for this system of equations 
must take into account three things: (i) the solutions 
must be geometrically regular at the origin; (ii) the so- 
lutions must yield an asymptotically flat spacetime; and 
(iii) the solutions must take into account the cosmological 
input for both the coupling a((p) as well as tp. 

Geometrical regularity at the origin means there is no 
conical singularity, i.e., the proper radius divided by the 
proper circumference should reduce to 2-7r at r = 0. This 
implies that g(0) = 1. Also, to maintain regularity in the 
field equations as r — > 0, we impose that d<&/dr\ r =Q = 
and d(p/dr\ r —o = 0. 

For a purely technical reason to set Mj = Me , we 
desire solutions that are asymptotically flat in both the 
Jordan and Einstein frames. That is, we want both g^ v 
and g^ u to reduce to the flat spacetime metric at spatial 
infinity. The implication of this is that the value of (poo = 
¥>(oo) must be such that 0(^00) = 0. This is guaranteed 
since there is one more rescaling that has no analog in 
GR. That is an invariance of the field equations if an 
arbitrary constant is added to the scalar-tensor coupling. 
If we simultaneously do the rescaling 



e x 



e c $ -> $ , e c N -> N , e c A 



A 



(2.29) 



on the variables defined by Eq. ( 2.21 ) and let a(ip) + c — ► 
a(tp), then the field equations remain unchanged. 

For the BD coupling, $ c = $(0) and <poo are the only 
freely specified field values. The value of N(0) is not 
specified freely, but rather is determined so that $(00) = 
0. The value of (p at the origin is not specified freely; it 
must be determined in such a way that the solution for 
tp goes to ifoo at spatial infinity. We will use the freedom 
to add an arbitrary constant to the BD coupling a(ip) so 
that all the solutions we consider have <p x = 0. 



C. The evolution equations 

We, here, assume the gravitational scalar field is time 
dependent Lp = ip(t, r) and use the rescaled bosonic field 
Analogous to eq.(3.6)-(3.10) in §, 
v and IT*: 



* as * = V4"7rG*V 

we introduce scalar field momenta n 



II V = j^d t (rip) = ^d t (rip), 



1 

T 
1 



II« = ift(rtf) = £$(rtt), 



(2.30) 
(2.31) 



where we set = N/g. The field equations become 
d t (np) = PIL V , (2.32) 



9 t n v = (d r /3)d r (np) + (3d r d r (rcp) - (d r 0)(r<p) 



1 



-Ngr 2ae 2a [-V CT *V CT * t + 2e 2a F(** f )], (2.33) 

$(rtt)=/m*, (2.34) 

d t n^ = (d r [3)d r (rt>) + j3d r d r (rt>) - (dr0)(r*)- 

r 

-2Ngr[e 2a ^ f J -aV g fV>]. (2.35) 



Note that \& and 11$ are complex variables, so ( 2.34 ) and 
( |2. 35 ) have two components. The momentum constraint 
and the G rr component of the Einstein equations become 

d tg = N[n v d r p + e 2a ^(uld r ^ + n»0 P *t)], 



N 



Nr 
T 



drN = --(.g 2 - 1) + ^ ( (dM 2 + n 2 - 



1 



We use the above set of equations ( |2.32 )-( p.37j ) for 
evolving the system and use the Hamiltonian constraint 
equation 



d r g g 2 - 1 _ n 2 



gr t 

„2a r n*nt 



(d r ipf 



+ ($.*)($.*+) + 2g 2 e 2a V {■*¥)] (2.38) 



to check the accuracy of our simulation. 



D. Numerical techniques 

1. Equilibrium configurations 

We use fourth-order Run g c-Ku tta algorithm to solve 
differential equations ( 2.23| )-(2.26). In order to find an 
equilibrium configuration, our system requires a two pa- 
rameter search to find a solution that satisfies the bound- 
ary conditions for both $00 and ipoo- Operationally, we 
choose a central value of the scalar field <& c first together 
with a guessed central value of the gravitational scalar 
field V'(O), and integrate out to large radii for different 
values of A(0). We then check if the resulting ipoo is 
close to our expected boundary value. 

The falling off behaviour for the gravitational scalar 
field ip is much slower than the matter scalar field <E>. 
Therefore, at the numerical boundary, say r = r enc i, we 
set the expected boundary value for (f(r e nd) &s f(fend) = 
</>oo + C/rendi where a constant C is given by C = 
. If the computed ip(r en d) is n °t 



„2 v d V (r) 
' end * dr 



(2.36) 



+e 2a [(d r V)(d r &) + II*l44 - 2g 2 e 2a V(^)} ] . (2.37) 



the expected value, then we change ip(0) and repeat the 
whole procedure. We set the tolerances to judge conver- 
gence in tp(r en d) as 5 x 10 -7 . More details are in mm. 



2. Evolutions 

We use the same code that was used in M for evo- 
lutions with modifications to incorporate BD theory. A 
polar slicing condition p6| for the lapse is hard wired 
into the code. This slicing is highly singularity avoiding 
and in the event of formation of an apparent horizon, 
the lapse rapidly collapses and the radial metric blows 
up crashing the code indicati ng im minent black hole for- 
mation. The lapse equation ( p. 37 ) is integrated once on 
every time slice using a sixth order i ntegr ation scheme. 
The Hamiltonian constraint equation (2.38) is monitored 
as an indicator of the accuracy of the simulation and is 
not solved during the evolution. A leapfrog evolution 
scheme is used as described in |]L0| . 

For the boundary conditions regularity dictates that 
the radial metric is equal to 1 at the origin. The boson 
field and the BD field are both specified at the origin. 
The boson field goes to zero at oo and the BD field goes to 
a constant which is fixed during the evolution. This con- 
stant does not enter into any of the evolution equations 
as all the terms in the set of equations are derivatives in 
the BD field. 

The inner boundary at the origin requires that the 
derivatives of all the metrics and fields vanish at this 
point. This is implemented by extending the range of 
r to include negative values. The metric compnents <?, 
N, the boson field and the BD field are required to be 
symmetric about r = 0. The code itself uses new field 
variables which are the original fields times r for both 
the boson fields as well as the BD field and the momen- 
tum variables II described in the previous section. Thus 



the code field variables as well as H v and n^, are anti- 
symmetric about r = 0. The bosonic and BD fields at 
the origin during the evolution are determined from the 
first derivatives of the code fields there. For the radial 
metric the value at the boundary is det ermined from the 
Hamiltonian constraint equation (2.38). 

In order to prevent reflections from the edge of the grid 
we use an outgoing wave condition at the outer edge of 
the grid (far away from the star) for the matter field W 
and the BD field ip. 

For the BD field, ip, we impose the wave behaves as 



1 



(p(r,t) =i Poo + -F{t--) 



(2.39) 



which is a natural out-going wave condition in spherically 
flat spacetime, and which is the proper assumption for 
the asymptotic region. Differentiating, one gets 



ldip 



dip 

dr 



((p-<Poo) 



0. 



(2.40) 



Note that this is the same technique used by Novak J25| , 
who studied a stellar collapse in scalar-tensor theory. 

For the boson field, 'J, an asymptotic solution of the 
form e~ kr e lQt /r to order - is assumed. This gives an 
outgoing boundary condition to this order of 



2 7„2 



I3 z k 



n 2 



N z m z e 



2 2a 



(2.41) 



This dispersion relation is nontrivial for a massive scalar 
field. There is no perfect algorithm to implement it. At 
the outermost grid point we require that 



dtdt* = -fidtdr* 



N 2 - 



^e z 



(2.42) 



where \& = r \&. The second term on the right is a finite 
mass correction to leading order. 

In addition to remove second order reflections we put 
a sponge for the matter field, (it reduces the momentum 
of the boson field artificially, and is irrelevant for the 
massless BD field) which is a potential term that is large 
for incoming waves (proportional to k + fi) and small for 
outgoing waves (proportional to k — CI). Therefore, we 
add a n additional term in the evolution equation for II^ 



e a N 



(n* + d r v) 



(2.43) 



for r e 



D < r < 



id where r en d is the r value of 



outer edge 



the outermost grid point and D is an adjustable param- 
eter representing the width of the sponge. D is typically 
chosen to be a few times the wavelength of the scalar 
radiation moving out. 

The code is tested with equilibrium configuration data 
(the zero perturbation case). We confirm that the evolu- 
tions leave the metric and BD fields perfectly static for 
several thousand time steps. We also checked the code 
with taking a large number of BD parameter lubd which 
would converge to the results in GR. 



III. EQUILIBRIUM SEQUENCES 

Before starting the dynamical study of boson stars, we 
describe briefly the equilibrium solutions in BD theory. 

The existence of ground state boson stars in BD the- 
ory are reported in ]17J for the case of ujbd = 6. The 
masses of the boson stars become smaller than the cor- 
responding configurations in GR. They also show the ef- 
fects of the interaction term, A. In [[L9|, the existence 
of higher node (excited) states are reported. They find 
that if the observational limit on ujbd is admitted then 
the obtained configurations are very similar to GR. They 
also describe the stability of the ground states config- 
urations using catastrophe theory. However, the whole 
discussion on the stability including the excited states is 
unconfirmed upto now. Also the dynamical behavior of 
this system are totally unknown. 
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FIG. 1. Stable configurations of a boson star in 

Brans-Dicke theory. Masses of ground state and excited state 
stars upto 3 nodes are shown as functions of the central mat- 
ter field "J> c . The solid lines and the dotted lines are the 



mass A / //(G»/m) 
spectively. 



and the particle number N p /{G*/m ), re- 



in Fig .pi, we show the sequences of equilibrium configu- 
ration of boson stars in BD theory both for ground states 
and excited states upto three nodes. These sequences are 



given by solving the set of equations described in § [I B , 



and we applied the BD parameter ujbd — 600 for plot 
ting this figure. The total mass of boson star M and 
its particle number N p (in unit m) versus central boson 
density $ c profile is plotted. We see that, on the ground 
and excited state sequences, the signature of the binding 
energy M — mN p goes from negative to positive according 
to the value of $ c . Configurations with positive binding 
energy are expected to be dispersive. 

In the case of GR, we know the star ground state con- 
figurations to the left of the central maximum are sta- 
ble. By 'stable', we mean that under perturbations they 
move to new configurations on the same branch (left of 
the maximum). We naturally suspect that the branch 
to the left of the maximum in the BD profile is also sta- 
ble. Therefore, in this and all further sections of this 
paper, we will refer to ground state boson star configu- 
rations to the left of the maximum mass in Fig.|lj as the 
5-branch configurations and others as [/-branch config- 
urations. From the catastrophe discussion |L9|, we also 
suspect that the [/-branch configurations are unstable. 

The boson star system exhibits anisotropies due to the 
presence of scalar fields. Here, "anisotropy" means the 
difference of the radial pressure from the tangential pres- 
sure in these configurations (\T£\ ^ |Tg|, where T" is 
the ener gy-m omentum tensor defined as the right-hand- 
side of (2.9). In the case of a scalar-tensor theory this 
anisotropy could change from GR due to the presence of 
the additional gravitational scalar field. In the case of 
BD theory, we find the anisotropy slightly but not signif- 
icantly lowers as compared to GR. In Fig.g, we show the 
fractional anisotropy 
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FIG. 2. Fractional anisotropy (3.1) is plotted for ground 
state configurations of central boson density $ c = 0.2 for 
GR and BD (ujbd = 600 and 60). The inset figure is the 
magnification of the range r — 7.0 to 7.4. 




FIG. 3. Fractional anisotropy (3.1) is plotted for 3 node 
(third excited) states in BD {ujbd = 600) with different cen- 
tral boson densities <& c = 0.05, 0.1, and 0.2. The plots end at 
the radius of each star (95% of the total mass). 



for a ground state boson star with <& c = 0.2 for GR and 
BD with ujbd = 600 and 60. In Fig.||, we show the 
fractional anisotropy for 3 nodes' (third excited) states 
of boson stars in BD ujbd = 600 with different $ c . At 
the nodes where the boson field derivative vanishes, all 
those fractional anisotropy are close to 2, which is the 
analytically expected number. This is because we can 
write the equilibrium fractional anisotropy as 



2[(d r tp) 2 +e 2a (d r $) 



(dr^) 2 



i Pr$) 



K$ 2 



2g 2 e 2a V(<5>)] 



(3.2) 



/„ = (t; - tDjt; 



(3.1) 



At the nodes <& = and as a result V($) also vanishes. 

The plots end at the radius of the boson star (we define 
the radius to be the 95% mass radius of these systems) . 
Consistent with Gleiser's results in the GR case [g8J, we 
have seen that the fractional anisotropy at the radius of 
the stars are almost the same number for all configura- 
tions. One of the consequences of these anisotropies is 
that we cannot apply the adiabatic perturbation method 



to discuss stability of this system, as discussed by Kaup 



IV. STABILITY OF GROUND STATE 
EQUILIBRIUM STATES 

In order to confirm whether the stars on the S-branch 
of the ground state in Figjlj are stable, we show our sim- 
ulation of the perturbed configurations of the S'-branch 
stars in this section. 

Our perturbations themselves simulate the annihila- 
tion (creation) of particles by enhancing (reducing) the 
boson field smoothly in some region of the star. A typical 
perturbation is put in for a region of the star from r\ to 
ri and is of the form 



f'2 



$perturbed{r) = <& original^) \ 1 + psin n (7T 



r 2 - r\ 

(4.1) 

where p is the perturbation size and is positive for particle 
creation and negative for annihilation, and n is an integer 
that can be varied. 

This perturbation is put in at the initial time step and 
the lapse equation and Hamiltonian equation are reinte- 
grated to give new metrics for this scalar field and original 
BD field. The momentum associated with the BD field 
is set to be zero on the initial time slice. 



A. Small perturbation of bosonic scalar field 

After we confirmed that our code retains its an equilib- 
rium profile for long time evolutions, we start our evolu- 
tions of an S'-branch equilibrium state with a tiny pertur- 
bation. Here, 'tiny' means that the difference of the mass 
between the original and new configuration is less than 
0.1%. We find that the system begins oscillating with 
a specific fundamental frequency, called a quasi-normal 
mode (QNM) frequency. Fig.[| shows the oscillations of 
the metric and the BD field for a boson star of central 
boson field density <I> C = 0.2 and with ujbd — 600. We 
see the BD field, as well as the metric function, start os- 
cillating at the same frequency. A Fourier transform of 
the metric and BD field data shows the overall frequency 
of the star to be - 0.03 in the nondimcnsional units. This 
is quite close to the case in GR. 

In GR, a QNM frequency increases as <& c become larger 
(radius become smaller) up to a point before starting 
to decrease rapidly to zero as it approaches the density 
corresponding to maximum mass signalling the onset of 
instability j8|. We found the same feature in BD. 

In addition, the system takes on the proper underlying 
frequency, which originates from the time dependence of 
the boson field \& ~ e lt . This frequency corresponds to 2-k 
period in t in our units, and appear in the metric and BD 




FIG. 4. Quasinormal mode oscillation of a stable boson 
star. The maximum value of the metric g rr , and the central 
Brans-Dicke field ip(r = 0) are plotted as functions of time. 
Both of them take on the QNM frequency of the star. The 
oscillation is virtually undamped for a long period of time. 



field oscillation with 7r period, from the structure of the 
equations. To show this feature we have enhanced the 
BD field oscillation (Fig.0) in Fig.pl. The time interval 
between the ticks on the horizontal axis is n. 




FIG. 5. The magnification of the second figure of FigM 
The time interval between the ticks on the horizontal axis is tt. 
We see the Brans-Dicke field has the underlying -k oscillation 
that it takes on. 



B. Large perturbation of bosonic scalar field 

Under large perturbations, a stable boson star in GR 
expands and contracts, losing its mass at each expansion. 
The oscillations damp out in time and the system finally 
settles down into a new configuration on the S-branch. 
These features are now observed in BD theory. 

We show the effects of a large perturbation on a sta- 
ble BD boson star described above. The example we 
present is the case of initial data with $ c = 0.2, mass 
M = 0.540G*/m after perturbation (about 13% lower in 
mass compared to unperturbed equilibrium configuration 
of M = 0.622G„/m). 

The maximum radial metric g rr and the central BD 
field as a function of time are shown in Fig.H and Figj^, 
respectively. In both figures, we plotted the case of 
<->Jbd = 600, 60 as well as GR. The increase of the max- 
imum g rr indicates the star is contracting, reaching its 
maximum value at the end of the contraction in a cycle. 
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FIG. 6. Finite perturbation of an 5*-branch boson star. 
Maximum metric g rr is plotted. The metric is damped in 
time as the star settles to a new configuration. 



-1.50 - 



/ v 


V 


\ 




to 


m =600 












1000 

lime 




1500 



V 



FIG. 7. The same model as FigM Central Brans-Dicke 
field tp(r — 0) is plotted. We see the oscillations damp out 
as the star settles, indicating again transition to a new stable 
boson star configuration. 



Then as the star expands, the maximum g rr decreases, 
reaching its minimum at the end of the expansion. These 
processes repeat themselves with the oscllations damping 
out in time as the star settles to a new stable configura- 
tion with maximum radial metric of smaller value than it 
started with (lower mass) . The lower value of the BD pa- 
rameter shows a phase shift in comparison to GR which 
might be suggestive of a different rate of approach to the 
final configuration. We see the same dynamical behavior 
in BD field (Fig.|^). It has the same oscillation frequency 
as that of the metric. The oscillations damp out in time 
and the BD field settles to a value closer to zero than it 
started at (lower final mass). 

The system loses mass through radiation during its 
evolution. A comparison of the mass as a function of 
time for BD case (u>bd = 600 and 60 as well as GR 
shows little difference, indicating that the radiation is 
mostly scalar field radiation and not scalar gravitational 
radiation. This is despite the BD field oscillating in the 
BD case and being zero in the GR case. 

The amount of mass radiated progressively decreases 



as can be seen in the luminosity profile (—dM/dt versus 
time) shown in Fig.[8|. Here again a comparison between 
GR and the two BD parameters is shown. Again we see 
the phase shift for the lower BD parameter (further from 
GR). The system finally settles to a new state of smaller 
mass. 

We have also conducted tests with perturbations that 
enhance the mass. For stable branch stars as in the case 
of GR S the star loses mass and settles to a new config- 
uration on the same branch. 
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FIG. 8. The same model as the previous figures. 'Lu- 
minosity' L = —dM/dt is plotted versus time. Clearly the 
radiation is decreasing in time. 



V. EVOLUTION OF [/-BRANCH STARS AND 
EXCITED STARS 

Boson stars on the [/-branch and excited states are in- 
herently unstable in GR [jnj]. Under perturbations that 
reduce the mass, boson stars on the [/-branch can mi- 
grate to the stable branch. In this section, we will see 
these features in BD theory too. They can be perturbed 
in a way so as to decrease their mass enough that they 
migrate to new configurations on the 5*-branch. On the 
other hand, in GR, if they do not lose mass and migrate, 
then the boson stars of configurations with M < N p m 
collapse to black holes. Stars with M > N p m are dis- 
persive and radiate out to infinity. We will see also these 
features in BD theory. 



A. Migration to Stable Configuration 

Our first dynamical example from [/-branch boson 
stars is a migration process. As in GR, we have also 
seen migrations of these stars to the stable branch when 
we remove enough scalar field smoothly from some re- 
gion of the star so as to decrease the mass by about 10%. 
In particular, we show the migration of a star of cen- 
tral boson field <& c = 0.35 with unperturbed mass 0.625 
G*/m. After perturbation, its mass is reduced to 0.558 



G*/m. In Fig|| we show the maximum value of metric 
g rr and central value of BD field ip(r — 0) versus time. 
The initial sharp drops in both lines occur as the star 
rapidly expands and moves to the stable branch. After 
that it oscillates and finds a new configuration to settle 
into. The damping of oscillations as it settles is clearly 
seen in the figure. We see also the BD field oscillations 
damping out as the star gets closer and closer to its final 
state. 
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FIG. 9. Migration of an unstable boson star to a stable 
configuration; central Brans-Dicke field ip{r = 0) and maxi- 
mum of g rr is plotted. There is a sharp initial drop in the 
radial metric as the star moves to the stable branch. The 
oscillations damp out in time as the star settles. 

The ratio of mass at time t to the initial mass for BD 
with parameter ljbd — 60, 600 and the GR case is shown 
in Fig.[t0| The flattening of the curve at later time is in- 
dicative of the star settling down to a new configuration. 
Although convergence towards GR with increasing u>bd 
is clearly indicated, there is no significant difference be- 
tween the three cases. The amount of the total mass ex- 
traction from the system is slightly supressed if we evolve 
in the BD theory. By the time of 7500 shown in the plot, 
we see that the mass of the star is about 0.045 G*/m, 
which corresponds to an equilibrium configuration with 
$ c = 0.06, while our central density $ c is about 0.061, 
meaning that the star is quite close to its final configu- 
ration. 



B. Transition to black hole 

Contrary to the previous example, if we add a small 
mass to [/-branch stars, we can see the formation of a 
black hole in its evolution. In Fig.|ll], we plotted an exam- 
ple of such evolution, which indicates formation of a black 
hole. The initial data of this plot is boson star of cen- 
tral boson field $ c = 0.35, which is the same value with 
the previous migration case with an unperturbed mass 
of 0.625G*/m). The system is perturbed very slightly so 
that the perturbed initial mass is 0.628G*/m, which is 
less than 0.5% greater than its initial mass. This is less 
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FIG. 10. Comparisons of total mass of the system M 
rescaled by its initial mass Mi n u during a migration process 
from (/-branch star. Three lines are plotted. Although the 
mass-loss is similar in the three cases the higher Brans-Dicke 
case clearly is closer to GR as it should be. 



than the maximum mass of a boson star of 0.632G»/to in 
the sequence of ujbd = 600. The sudden collapse of the 
lapse function is indicative of the imminent formation of 
an apparent horizon due to the polar-slicing condition in 
our code [£6| . In addition to this the radial metric starts 
to blow up and the code is no longer capable of dealing 
with the sharp gradients and crashes. As an indicator 
of the suddenness of the process, we see that in the con- 
figuration shown the lapse has fallen to a value of about 
0.003 by a time of 60 after being at 0.230 at a time of 
55 and 0.5 at a time of 50 (the latter two points are not 
shown in the plot). The lapse value at the edge of the 
grid is greater than one because it has been scaled by the 
underlying frequency of the system f2. 

There is almost no loss in mass in this system and the 
time of collapse is quite similar to GR. In the GR case, 
we confirmed the formation of black hole [p9f very shortly 
after this point («3M where M is the mass of the sys- 
tem) by switching this data into 3-dimensional code and 
evolving the system. Therefore we expect almost the 
same behavior in the BD case. 

Note that we have seen this behavior for denser boson 
stars with lower masses and the same order of perturba- 
tion. This is true as long as they have M < Nm. 

The black hole formation in the BD theory is investi- 
gated by Scheel, Shapiro and Teukolsky for dust collapse 
case |23|. They found that the dynamical behavior of 
horizons (both event horizon and apparent horizon) are 
quite different in the physical Brans-Dicke frame, but the 
same as GR in the Einstein frame. Their results there- 
fore also support our discussion of the formation of black 
hole in the BD theory. 

We can not show the emitted scalar waveform because 
we are using polar slicing condition and can not pro- 
ceed with the evolution after the collapse of lapse. We 
are planning to evolve this system after the formation 
of black hole in 3 dimensional code. The results will be 



reported in the future. 
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FIG. 11. Dynamical transition from [/-branch star to a 
black hole. The metric poo is plotted. The collapse of the 
lapse function is indicative of imminent black hole formation. 



C. Transition to ground state from excited state 

Excited states of boson stars in general are not stable 
in GR case. They form black holes if they cannot lose 
enough mass to go to the ground state. We confirm the 
same features for the BD case. In Fig.O, we plot the 
metric g rr of dynamical transition from an excited state 
with 3 nodes to a ground state boson star configuration. 
This star had an unperturbed mass of 3.249 G*/m cor- 
responding to a central density $ c = 0.01. After pertur- 
bation its mass was reduced to 0.919 G*/m. The initial 
configuration has four metric maxima and the final has 
one showing the transition. After it goes to the ground 
state, it oscillates and compactifies to form a new config- 
uration. We show the oscillations of the star from a time 
of 27300 to 28400 in Fig ||. The 95% mass radius at this 
stage is about 100 and the star has still to continue its 
evolution for a while longer. The amount of mass loss 
in this process is quite similar to the migration case (cf. 
Fig.|Ty) . We also found that the difference of the theory 
is little. 

Any configuration of an excited state star with a mass 
less than the maximum mass of a ground state star is ex- 
pected to go to the ground state. However since higher 
excited state S branch configurations have progressively 
greater masses for the same central density ( Mass of n- 
nodes > n — 1 > ■ • • > 1 > ground state boson stars) 
the configurations with masses less than the maximum 
ground state mass get more and more dilute. We have 
seen the tendency of these stars to migrate to the ground 
state in our tests. In ]l(J, such cases have been described 
for the GR case. Since these dilute configurations cor- 
respond to large stars with low oscillation frequencies, 
to evolve them until they go to the ground state is nu- 
merically very costly as they take a very long time to do 
so. Configurations with slightly larger masses may also 



be able to lose enough mass during their evolution and 
transit to the ground state. An exact parameter search 
though would be very time consuming. We have rather 
chosen to just exemplify the possibility of transitions with 
our example. We have taken a denser configuration to re- 
duce the time scale of evolutions and then perturbed it so 
much to ensure that no black hole forms. In fact we can 
just think of the system as a perturbed distribution of 
bosonic fields that has 3 nodes and which transits to the 
ground state rather than compare the perturbed system 
to the original unperturbed one. 




FIG. 12. Dynamical transition from an excited state to 
a ground state boson star configuration. The metric g rr is 
plotted. The initial four peaks indicative of a 3-node star 
cascades to the ground state after a long time evolution. 
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FIG. 13. Dynamical transition from an excited state to 
a ground state boson star configuration. The metric g rr is 
plotted at later times to show its oscillations after the star 
reaches a ground state. 



VI. FORMATION OF BOSON STARS IN 
BRANS-DICKE THEORY 

In the previous sections, we analyzed boson stars in 
BD theory, starting with those equilibrium or perturbed 
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equilibrium configurations. However, we have not dis- 
cussed whether or how such an equilibrium configuration 
actually forms in BD theory. In this section, we answer 
this question by demonstrating the formation of boson 
star in BD theory. The formation of boson stars in GR 
has been discussed by Seidel and Suen [p0| . 

We start our evolution with the initial data that has a 
Gaussian packet in the bosonic field $, which represents 
a local accumulation of matter field: 



$ = acxp(— bx 2 ) 



(6.1) 



where a and b are free parameters. We set the BD field <f 
to be flat at the initial stage, so as to see if local inhomo- 
geneity of the matter will form a b oson star in BD theory. 
We integrate the lapse equation ( [2.36J ) and Hamiltonian 
constraint equation (2.38) to provide metric functons on 
the initial slice . We, then, evolve the set of dynamical 
equations (|2.32|)- (|2.35| ). 



boson oscillations corresponds to masses between 0.342 
and 0.410 respectively while the BD field oscillations give 
a mass between .355 and .405 respectively. Given that 
the mass at this stage is 0.384 (consistent with the above) 
we expect that the final mass will be between 0.355 and 
0. 
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FIG. 14. An example of formation of boson star in 
Brans-Dicke theory. Snapshots of fBD (r) are plotted for ini- 
tial stage of evolution. 

We find that, with particular parameters a and b, this 
system actually forms a stable, equilibrium configuration, 
which might be recognized as the formation of a boson 
star. As a demonstration, we here show an evolution 
with parameters a — 0.1 and b = 0.025. The BD param- 
eter ubd is taken to be 600. In FigJlJ, we show the BD 
field if as a function of radial coordinate at various ear- 
lier times of evolution. We see that the BD field becomes 
negative quickly and begins oscillating around a partic- 
ular value. Actually the reader will find the BD field 
ip is jumping around r = 40 — 50 at t = 50. The long 
time evolution of BD field is shown in Fig.flq, in which 
we show the BD scalar field ipso at the center for earlier 
times, middle times and later times. We can see the field 
settling down to a periodic oscillation in final phase, like 
in the migration and transition cases in the previous sec- 
tion. The initial mass of this configuration in units of G* 
is 0.39 and the final(at the end of our simulation) about 
0.384. At this stage the magnitude of the central boson 
field is oscillating between 0.032 and 0.048. The BD field 
oscillates between —0.00126 to —0.00166. This range of 



FIG. 15. An example of formation of boson star in 
Brans-Dicke theory. Dynamical behavior of the Brans-Dicke 
scalar field ipbd(x — 0) is plotted for three evolution regions: 
earlier time, middle time and later time. We can see the field 
settling down to an equilibrium configuration (periodic oscil- 
lation) . 

We show the luminosity L(= —dM/dt) versus time 
t curve in Fig.Qq within the same periods as in Figjlq. 
In the early stage, we see one pulse is emitted from the 
system. This is related to the outgoing pulse from our 
initial boson field setting. After this initial pulse the 
system slowly takes on the characteristic w oscillation of 
l^l 2 of the star as the star begins forming. After that, 
we see L begins damped oscillations versus t. (We cut 
out the initial large amplitude luminosity around t — 
200). The system's evolution is followed for a long time, 
however, the accuracy of the calculation is quite good 
with the hamiltonian constraint satisfied to order 10~ 7 
or better. 

We also note that certain parameters a and b (men- 
tioned at the beginning of this section) will result in for- 
mation but others will not. If we choose large amplitude 
a and small &, then the initial configuration has too large 
a mass and is not dispersive enough resulting in black 
hole formation during evolution. In the opposite limit 
if we have a very narrow localized wave packet it has a 
tendency to be dispersive (wave equation). So if b is too 
large no boson star forms. Intermediate between these 
two are configurations that form stable stars. For ex- 
ample, if a = 0.1, the b — 0.01 a black hole is formed, 
b = 0.025 results in boson star formation as we have 
shown, and b — 0.035 become flat space at the end of the 
evolution. On the other hand for b = 0.01 and a = 0.05 
a boson star forms. 
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FIG. 16. Similar to FigJlS, the emitted luminosity 

L = —dM/dt is plotted. In the first evolution stage, the 
luminosity data takes on the underlying boson field square 
oscillation, after having emitted one scalar pulse, related to 
the initial field configuration: Gaussian pulse. The amount 
of mass loss decreases in time as the formed star settles. 



VII. CONCLUDING REMARKS 

We studied dynamical features of boson stars in the 
Jordan-Brans-Dicke (BD) theory of gravity By evolv- 
ing the system numerically, we discussed the stability of 
those ground state S'-branch equilibrium configuration, 
black hole formation and migration from [/-branch solu- 
tion, and transition process from excited state to ground 
state, together with those scalar wave emissions (mass- 
loss from the system). 

We showed that the basic features are the same as the 
general relativity (GR) cases |1(J. Since we choose BD 
parameter ljbd = 600 for most simulations as we are in- 
terested in those observationable differences with GR, the 
signatures of boson stars such as oscillation frequency of 
the final S'-branch configurations and the emitted waves' 
luminosity are quite identical with GR. The scalar wave 
emissions are slightly supressed than GR. These results 
indicate to us that the scalar modes of gravity in the bo- 
son star system is already dominated by the boson scalar 
field, thus additional gravitational scalar field does not 
change the dynamical behaviour in BD theory of gravity. 

We also demonstrated the formation of a boson star 
from a Gaussian packet of bosonic field and flat BD field. 
This success suggests to us that the boson star is real- 
izable object even in the BD theory, opens windows to 
study them in astrophysical roles including similar non- 
topological solitonic objects. 

By showing that the [/-branch of ground state equilib- 
rium is unstable, we have confirmed the stability discus- 
sions by Comer and Shinkai for BD ground state using 
catastrophe theory. Our next problem is to check their 
prediction in the scalar-tensor theory, especially Damour- 
Nordtvedt attractor model |2(j] , with which they predict 
that the boson stars in the early universe will not be 
formed. We are also now studying the existence of "oscil- 



latons" pi), a self-gravitating solitonic object made from 
a real scalar field, in the scalar-tensor theory. These re- 
sults will be reported elsewhere. 
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